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50103-422 

REAL TIME STATISTICAL COMPUTATION IN 
EMBEDDED SYSTEMS 

Related Application 

[0001] This application claims the benefit of U.S. Provisional Application No. 

60/283,779 entitled "Real Time Statistical Computation in Embedded Systems" filed on April 
12, 2001, the disclosure of which is entirely incorporated herein by reference. 



5 Field of Invention 

[0002] The present subject matter relates to an efficient technique for computing standard 

p deviation (sigma) and variance for a variable, such as a variable parameter measured as part of a 
)~ closed-loop control operation or a test procedure, and to circuitry and/or programs and 
□ processors for implementing that technique. 
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10 Background 

□ 

fly [0003] There are many applications today that involve monitoring of a physical process 

j'pi 

jjjjj and computing useful statistical information based on sample data obtained from the process 
h3 monitoring. In a manufacturing process, this may entail sampling a process variable or a 

ru 

performance variable of the manufactured process, for example as each unit produced passes 
15 through a particular phase of the production process. In an electrical or mechanical device, this 
may entail sampling a detectable parameter during device operation, for example as part of a test 
run immediately following manufacture or during a later re-test. In a disk drive, for example, the 
drive processor may detect and process a temperature value or a "seek time" required to find an 
addressed location on the disk. In a tape drive example, the drive processor may monitor a 
20 position error signal (PES), the Tracking Servo Error, the Tape Speed Error, under run margin, 
under run time or ramp distances. 

[0004] Running control functions often simply determine the difference between a 

current value of such a measurement and a target or threshold value and adjust the process in 



question based on the result of that determination. However, there are many testing and 
manufacturing control techniques that require more complex statistical information. 
[0005] In many modern processes and devices, it is necessary to control a manufacturing 

process or other operational process based on a measure of performance over some period of 
time. For this purpose, the circuitry or processor must process some number of samples and 
compute a useful measure of performance from the sample population. A large number of such 
applications call for statistical calculation of sigma (a), which is the standard deviation; and 
other processes utilize the variance, which corresponds to cr 2 . Other performance measures are 
calculated from a or a 2 , such as the capability analysis index (Cpk) of a process, which shows 
whether process data points are within specified limits. The process capability index, Cpk, is a 
standard measure of process capability over an extended period of time for a process exhibiting 
statistical control. Cpk is considered to be a reliable indicator of process performance, taking into 
account process variation and deviation from nominal. A "design of experiment" (DOE) process 
changes a process variable, and the a, the a 2 , or some related performance measure is calculated 
during a subsequent process run to determine the process performance. 

[0006] For example, US Patent No. 5,956,251 to Atkinson et al. discloses methods for 

meeting end item/assembly tolerance criteria for large flexible parts, which includes calculation 
of values of Cpk based on the standard deviation (a) of the sample population. US Patent No. 
6,269,326 to Lejeune teaches processing the population of samples for each of a plurality of tests 
run at different measurement dates, to compute an average or mean of the measurements and a 
standard deviation of the measurements (sigma) for each test date. Then, for each set of samples 
and for each test date, the methodology computes a criterion of appreciation (CP), defined as the 
ratio between a difference of limits and the standard deviation. 

[0007] In current processing and control techniques, the calculation of a or a 2 requires 

considerable time to accumulate sample data and excessive data storage for the samples, until the 
a or a 2 value(s) can be calculated using an entire sample data population. Even after collecting 
the entire population, the individual processing of each sample as part of the computation takes 
considerable time and processor power. More specifically, calculation of a or a , requires 
accumulating the entire "population" or sample set and then calculating the mean (sum divided 
by n, the number of samples). The mean is subtracted from each individual sample data point, 
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and each result is squared. The variance (a ) equals the sum of these squares divided by n-1; and 
of course, the standard deviation (ct) equals the square root of the variance. 
[0008] For an n-sample size, the collection of data can consume a sizable memory and 

the processing thereof requires considerable time after the data acquisition is completed. As a 
result, real-time computation of the a or a 2 value or of any of the performance measures that 
may be derived therefrom is not practical. 

[0009] Hence, a need exists for a particularly efficient technique for computing standard 

the deviation (sigma) and/or the variance of a variable in real time, for example, for use in a 
process controller or in a performance measurement application. A related need exists for a 
processing circuit and/or a programmed general purpose processor for implementing the efficient 
technique for computing these statistical measures of the variable. 

Summary 

[0010] The inventive concepts meet the above noted needs in the art as they relate to 

techniques and devices for efficiently calculating standard deviation or variance or the like, for a 
measured process variable or parameter. For example, the concepts alleviate the need to store all 
samples in the test population until the end of the test run, substantially reducing the memory 
requirements. The concepts also provide a simpler calculation at the end of the test run, reducing 
the processing power and time required. These improvements make it possible to implement the 
calculation of standard deviation or variance for a measured variable or the like as an imbedded 
function of a processor or other circuit within a manufactured electromechanical device, such as 
a tape drive or disk drive. 

[001 1] Hence, one aspect relates to a device for computing a statistical value related to a 

measured process parameter during ongoing operation of the process. A sampling circuit, 
responsive to a signal representing the measured process parameter, provides a predetermined 
number of samples of the measured process parameter in sequence, during operation of the 
process. The device includes means for generating a representation of the variance of the 
measured process parameter and/or the standard deviation of the measured process parameter in 
real-time. The device does not store all of the samples of the measured process parameter until 
completion of the sampling. In a preferred embodiment, the means generates the representation 
based on summation of the predetermined number of the samples of the measured process 
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parameter and summation of squares of the predetermined number of the samples of the 
measured process parameter, while actual sample values are effectively discarded. 
[0012] A method embodiment in accord with this concept generates a statistical measure 

of performance, such as variance or standard deviation, from a measured process variable during 
5 ongoing operation of a process. This method entails measuring the variable to generate a signal 
and taking a predetermined number (n) of samples from the signal, during the ongoing operation 
of the process. Concurrent with the sampling, a running sum of the n samples is accumulated, 
and a running sum of squares of the n samples is accumulated. Accumulation of these two sums 
does not require storing all n samples until the end of the method. Typically, the individual 
10 samples can be discarded. At the end of the sampling, the method involves processing a final 
value of the sum of the n samples and a final value of the sum of the squares of the n samples, to 
produce the statistical measure of performance. 

[0013] An embodiment of an apparatus for implementing this technique includes a 

J*} sampler responsive to the signal representing the measured process variable, for sampling the 

i.iy 

□ 15 signal during ongoing operation of the process to generate a predetermined number (n) of 

\'f\ 

in samples. The apparatus includes two computation modules, which may be implemented in 
M " * circuitry, but preferably are implemented in code or microcode for execution by a programmable 

» 

□ digital processor. An interim computation module accumulates the sum of the n samples, and it 

HI 

accumulates a sum of squares of the n samples, during the ongoing operation of the process, 



;;Jj 20 without retaining all of the n samples. A one time computation module, coupled to the interim 
f(J computation module, computes the statistical measure of performance in response to the final 
sum of the n samples and the final sum of squares of the n samples. 

[0014] The preferred embodiment of the interim computation module comprises two 

accumulator loops and a multiplier. A first accumulator loop receives the n samples and 

25 accumulates the sum of the n samples during ongoing operation of the process. The multiplier 
calculates the square of each of the n samples as received from the sampler, and the second 
accumulator loop accumulates the sum of the squares of the n samples from the multiplier. 
[0015] The preferred embodiment of the one time computation module comprises means 

for receiving the sum of the n samples and the sum of squares of the n samples and in response 

30 calculating: 
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[0016] 2 1 f ■ 2 2 > 



[0017] where: 



[0018] a 2 is variance, 

[0019] ji is mean of the n samples of the process variable (sum of sample divided by n); 

5 and 

[0020] Xj is sample for ith sampling interval, in range from 1 to n. 

[0021] In the preferred embodiment, the one time computation module also outputs the 

mean (n) and computes the square root of the variance to provide an output of the standard 
deviation. 

10 [0022] A particular advantage of the technique supplemented by the embodiments is that 

u 

□ it can be implemented as an embedded function of any manufactured device that contains an 
appropriate electronic circuit or microprocessor. The manufactured device can include discrete 

O logic circuitry or preferably one or more modules of program code for the processor, to enable 

i!|f= 

i.p the device to automatically calculate the variance and/or the standard deviation and possibly 

rif? 

'"" 15 other statistical measures of performance, for example as part of a test run immediately 

ii 

□ following manufacture or during a later re-test. In a disk drive, for example, the drive processor 
j'pj may detect and process a temperature value or a "seek time" required to find an addressed 

location on the disk. In a tape drive, the drive processor may monitor a position error signal 
y (PES), the Tracking Servo Error, the Tape Speed Error, under run margin, under run time or 
20 ramp distances. In each of these examples, after accumulating the sum and the sum of the 
squares from the desired number of samples, the imbedded circuit or program functions to 
perform the one-time calculation(s) and output one or more of the statistical measures of 
performance. 

[0023] Additional objects, advantages and novel features of the embodiments will be set 

25 forth in part in the description which follows, and in part will become apparent to those skilled in 
the art upon examination of the following and the accompanying drawings or may be learned by 
practice of the invention. The objects and advantages of the inventive concepts may be realized 
and attained by means of the instrumentalities and combinations particularly pointed out in the 
appended claims. 
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Brief Description of the drawings 

[0024] The drawing figures depict preferred embodiments by way of example, not by 

way of limitations. In the figures, like reference numerals refer to the same or similar elements. 
[0025] Fig. 1 is a flow-chart, essentially in block diagram form, showing the high level 

5 functional elements or steps of an exemplary embodiment. 

[0026] Fig. 2 is a flow-chart, essentially in block diagram form, showing the processing 

elements or steps of the interim (Int) computation module in the embodiment of Fig. 1. 
[0027] Fig. 3 is a flow-chart, essentially in block diagram form, showing the processing 

elements or steps of the one time computation module in the embodiment of Fig. 1. 
10 [0028] Fig. 4 is a simplified functional block diagram of a data drive, for example a disk 

drive or a tape drive, having the exemplary embodiment of Figs. 1-3 embedded therein. 

Detailed Description of the Preferred Embodiments 



q [0029] For various analyses, it is useful to compute the mean of the data, since the mean 

■: s [ is the simplest computation of a measure of performance that does not require acquisition of the 
CP 15 full set of sample data prior to the calculation. The mean is the sum of all data divided by the 

IS 

total sample size. Hence, a simple algorithm samples the data and accumulates both sample size 
j t jj (count) and the sum of the sample data, and this algorithm can continually calculate the mean by 

[fl simply dividing current value of the sum by the current vale of the sample size (count). 

i-3 

jipi However, the final mean is calculated only when the entire population of sample data is finally 

20 accumulated and processed. 

[0030] Variance is the square of the standard deviation. Computation of the variance 

requires both the square of the individual sample values and the mean of the process. Since we 
cannot get the mean until the total sample is collected, the online computation of variance as we 
acquire samples becomes a non-trivial task, without collecting and later processing the entire set 

25 of sample data. Of course, it is desirable to compute both the means and the variance without 
waiting to complete the collection of the entire set of sample data, since such collection requires 
a large buffer side for the process and requires delay during accumulation of all of the samples. 
The inventive approach therefore calculates both the mean and the variance of a given variable, 
but without having to collect the entire N size set of samples of the data. 
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[0031] Instead of storing the entire set of samples and performing all the calculations of a 

desired statistical measure of performance from the sample set, the mean (based on a running 
summation of the samples) and a running summation of squares of the samples are used. Hence, 
the embodiments disclosed herein relate to a technique for efficient calculation of one or more 
5 statistical measures of performance based on samples of a monitored parameter or variable of an 
ongoing process. The inventive technique involves accumulating running totals of: (1) the 
samples of the variable and (2) the squares of the samples of the variable. However, it is not 
necessary to store either the individual samples or the individual squared data points. Essentially 
both the samples and the squares are discarded. When all the sample values have been taken, 
10 and the two summations are complete, simple one-time calculations based on those two 
summations can quickly and accurately provide the variance and/or the standard deviation. The 
preferred embodiment also provides a calculation of the mean value. The one time computation 
q can provide any or all of the statistical values without using the individual samples or their 
J jj individual square values. Other performance measures may be quickly derived from one or more 

□ 1 5 of these three values. 

lip 

i'g [0032] As discussed later, the inventive embodiment relies on a new formulation for 

i,M calculating variance and/or standard deviation as measure(s) of performance. To understand 

ii 

□ how this is possible, it may be helpful first to consider the following mathematical derivation. 

ru 

py [0033] It is well settled that the standard deviation a is defined by the following formula: 



!:Q 20 




[0034] o- = iU (1) 



[0035] where n is the number of samples of the variable data x. The value Xj is the 

sample value of the variable x taken at the ith point in time or sample interval. The expression 
shown above at (1) is the normal formulation used to calculate the standard deviation. 
25 [0036] The mean fa, which is the sum of the sample values divided by the number of 

samples, can be represented by the following formula. 
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[0037] A = (2) 

n 

[0038] The variance is essentially the square of the standard deviation (a). Hence, 

applying expression (1), the variance can be expressed as shown below. 

[0039] w-ltT (3) 

[0040] The expression (3) for the variance involves a summation of square terms, where 

each the terms is the difference between the ith sample (xj) value and the mean (u). The actual 
multiplication of the squaring function can be expressed as (xi - u.) * (x; - When this 
multiplication is carried out, within the expression (3), the expression for the variance becomes. 



[0041] 



a 2 =- 



n-\ 



2>,. 2 -2//**,+// 2 ) 



(4) 



[0042] It is then possible to distribute the summation function to all of the terms within 

the innermost expression. 



[0043] 



a 2 = 



f n 



n-\ 



n n \ 

2 



Vi=l i=l M J 



(5) 



[0044] However, since the mean (jli) is a constant, the summation of the mean over n 

sample intervals becomes the product of the mean times the number n of the samples, which 
stated as an equation, reads: 

[0045] 2> 2 =»A 2 - (6) 

1=1 

[0046] Applying (6), the expression (5) for the variance becomes: 
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[0047] 



a 2 = 



\ 1=1 i=i J 



(7) 



[0048] Recall that from expression (2) the mean (n) equals the sum of the sample data 

values Xi divided by the number n of samples. From that, we know that we can: 



[0049] 



n n 

replace^ x { with J] x ( = nju 



1=1 



5 [0050] When we make this replacement in the middle term of the expression (7), then the 

expression for the variance becomes: 



[0051] 



a 2 =- 



1 



f n 



n-\ 



Y^x 2 -iM-n-ju + npi 1 



(8) 



sD 
® 

!i 

hoi 



[0052] Then, the expression can be simplified further by performing the multiplication of 

the |i terms, so the expression becomes: 



10 [0053] 



a 2 = 



1 



( n 



n-\ 



^x 2 -2n M 2 +n M 2 



\i=\ 



(9) 



[0054] Then, the expression can be simplified by combining like terms, so the expression 

for the variance becomes the result: 



[0055] 



a 2 = 



n~\ 



(10) 



15 [0056] Unlike the expression (1), the resulting expression (10) can be calculated from 

data accumulated in real-time, during processing, without storing all of the original sample data. 
In this expression (10), the n term is a constant, that is to say the desired number of samples. For 
calculation of the mean (u), it is necessary to accumulate the total of the sample values, however, 
this can be done on a running bases during operation of the processes. For the square summation 
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term, it is possible to square each sample value x; when received and add the result to a running 
total in a summation or accumulation register, during the process operation. When all the 
samples have been processed, and the totals of these two values (x and x 2 ) have been 
accumulated, it is only necessary to calculate the mean from expression (2) and perform a one 
time calculation of the result shown in expression (10) to obtain the variance. A square root 
calculation will then yield the standard deviation (cr). Various desired measures of performance, 
such as Cpk, can be derived from the mean, the variance and/or the standard deviation. 
[0057] Reference now is made in detail to the presently preferred embodiments for 

implementing this technique, examples of which are illustrated in the accompanying drawings 
and discussed below. Fig. 1 is a block diagram of an implementation of the technique outlined 
above. Figs. 2 and 3 provide similar diagrams for two of the elements of the embodiment of Fig. 
1. The illustrated implementation may utilize discrete logic components, essentially as shown. 
With the possible exception of the physical sensor and associated sampler circuit, the functional 
blocks preferably are implemented in code executed by a processor. For example, in a tape drive 
or a disk drive, most of the elements shown in Fig. 1 may be implemented in microcode run by 
the drive processor. However, for ease of discussion and understanding, the elements shown will 
be referred to as elements or modules. 

[0058] Fig. 1 shows the elements of a circuit 10 for computing the mean (fi), the variance 

(a 2 ) and the standard deviation (cr) of a set of samples of a measured process variable. In 
practice, the circuit 10 will include a sensor 11 or other element for detecting the variable 
parameter. In a disk drive, for example, the sensor may be a temperature sensor. If the sensor 1 1 
provides an analog signal, then the circuit 10 will include a sampler 13. The sampler may be as 
simple as an analog to digital converter (ADC), for example if the sensor 11 is an analog 
temperature sensor. In other applications, the sampler 13 may process one or more measured 
signals to generate samples of a desired parameter, such as a position error signal (PES), the 
Tracking Servo Error, the Tape Speed Error, under run margin, under run time or ramp distances. 
The sensor 11 typically is a hardware element, and the sampler 13 may include some hardware, 
although in some embodiments some or all of the functions of the sampler may be performed by 
execution of program code by an appropriate digital processor. 

[0059] A pulse generator 15 provides clock signals at the desired sampling rate, for 

example, to the sampler 13 to control the timing and rate of the analog to digital conversion. The 
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pulse generator may be a separate clock circuit but preferably is implemented as an element of 
the processor controlled by appropriate program code. 

[0060] The sampler 13 provides a sequence of samples of the sensed variable signal to 

the input (Inl) of an interim computation module 17, at the rate controlled by the sampling clock 
signal from the pulse generator 15. The module 17 also receives the clock signal from the pulse 
generator 15. The module 17 is a digital circuit or a process routine in a programmed digital 
processor for calculating two specific sums in response to the input samples and the sampling 
clock. In the embodiment, the interim computation module 17 accumulates the two summation 
values and provides those values at output ports Outl and Out2. The first output (Outl) is the 
sum of the input samples. The second output (Out2) is the sum of the squares of the input 
samples. The interim computation module 17 supplies these two output values to appropriate 
inputs of a one time computation module 19. 

[0061] The circuit 10 also includes a counter 21, coupled to the pulse generator, for 

counting a predetermined number of samples or sampling intervals. Essentially, the counter 
provides a control signal to the one time computation module 19 when the circuit has received 
and processed the desired number n of the samples of the variable parameter, to trigger the one- 
time computation by that module. The module 19 also receives the number n, as a constant 
value indicating the number of samples (on Inl), for example from the circuit 23. 
[0062] Hence, the one time computation module 19 receives the number n of sample 

values on its first input (Inl), the summation of the n samples on its second input (In2), and the 
summation of the squares of the samples on its third input (In3). When triggered by the counting 
of the n samples by the signal from the counter 21, the one time computation module 19 
computes the mean as the value of the sum of the samples divided by the number n of the 
samples, and outputs that value on its first output line (Outl). Concurrently, the one time 
computation module 19 computes the variance (a 2 ) based on expression (10) above, and it 
computes the standard deviation (a) as the square root of the computed variance. 
[0063] Hence, the illustrated embodiment of the one time computation module 19 outputs 

the mean (u), the variance (a 2 ) and the standard deviation (ct) immediately following the taking 
and processing of the last of the n samples. The module 19 could provide these computed 
statistical values to one or more displays. In an embodiment embedded in a disk or tape drive or 
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other manufactured device with a processor, the processor incorporating the module 19 would 
output the three statistical values via one or more appropriate data ports. 

[0064] To insure a full understanding of the embodiment, it may be helpful to consider 

examples of the interim computation module 17 and the one time computation module 19. Fig. 2 
shows a possible embodiment of the interim computation module 17, and Fig. 3 shows a possible 
embodiment of the one time computation module 19. For convenience, the elements for 
triggering the operations of these modules, in response to the clock or the counter output, have 
been omitted, although various means for implementing such timing/control functions should be 
familiar to those skilled in the art. 

[0065] As shown in Fig. 2, the samples received on the input (Inl) are periodically 

applied to a first input of an adder 25 of the module 17. The output of the adder 25 is supplied to 
the input of an accumulator register 27, which provides a delay (1/z) equal to one sampling 
interval. The register output is fed back to the second input of the adder 25, so that each new 
sample is added to the output register value (which corresponds to the sum produced in the 
previous sampling interval). In this manner, the adder 25 and the first register 27 form a first 
accumulator loop, for processing the sample values. Over a number n of sampling intervals, the 
register 27 accumulates the total sum of the set of n samples of the input variable and provides 
this total on the first output line (Outl) of the interim computation module 17. 
[0066] The samples received on the input (Inl) also are periodically applied in parallel to 

both inputs of a multiplier circuit 29 or other two-input product device. During any given 
sampling period, the multiplier circuit 29 produces the square of the current sample. 
[0067] During each interval, the multiplier circuit 29 supplies the square of the current 

sample to a first input of an adder 31. As the samples are received in sequence, the multiplier 
circuit 29 sequentially supplies the squares of the samples to this input of the adder 31. The 
output of the adder 31 is supplied to the input of a second accumulator register 33, which 
provides a delay (1/z) equal to one sampling interval. The register output is fed back to the 
second input of the adder 31, so that each new square of a sample is added to the output register 
value (which corresponds to the sum of squares produced in the previous sampling interval). In 
this manner, adder 31 and the register 33 form a second accumulator loop, which processes the 
squares of the samples. Over a number n of sampling intervals, the register 33 accumulates the 
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total sum of the set of n squares of the samples of the input variable, and the register 33 provides 
this second total on the second output line (Out2) of the interim computation module 17. 
[0068] The one time computation module 19 (Fig. 3) uses the sum of the samples, which 

it receives on its second input (In2) to calculate the mean, for example in accord with expression 
5 (2). In accord with the embodiment, the one time computation module 19 uses the sum of the 
samples from its second input (In2) as well as the sum of the squares received on its third input 
(In3) to calculate the variance, for example in accord with expression (10). In the embodiment, 
the module 19 computes the standard deviation as the square root of the computed variance. 
[0069] The computation module 19 is triggered to perform its computations at the end of 

10 n sample intervals, that is to say when a complete set of n samples have been taken by the 
sampler 13 and have been processed by the module 17. As shown in Fig. 3, value n (from input 
Inl) is applied to an inverter 35, which provides a digital value representing 1/n. The value n 
(from input Inl) also is applied to a circuit 37, which provides a digital value representing l/(n- 
1). The circuit 37, for example, may comprise a subtraction circuit for calculating 1-n and an 

1 5 inverter for calculating the inverse of the subtraction result. 

m 

!.p [0070] A multiplier circuit 39 or other two-input product device receives the value of 1/n 

m 



a 
o 



and multiplies that value by the total sum of the sample values received on the second input (In2) 
□ of the module 19. The product output by the circuit 39 is the mean (ja), which is supplied to the 

ru 

m first output line (Outl) of computation module 19. 

p 20 [0071] The multiplier circuit 39 also supplies the mean (ji) in parallel to both inputs of a 

MJ multiplier circuit 41 or other two-input product device. The multiplier circuit 41 therefore 

produces a product value equal to the square of the mean value or equal to u 2 . 

[0072] The multiplier circuit 41 supplies the value \x to one input of another multiplier 

circuit 43 or other two -input product device. On its other input, the multiplier circuit 43 receives 
25 the constant value n (from input Inl). Hence, the product generated by the multiplier circuit 43 

equals njo. . 

[0073] On its positive input, a subtractor 45 receives the sum of the squares of the set of 

samples from the interim computation module 17, via the third input (In3) of the module 19. 
The subtractor 45 receives the value of nu 2 from the multiplier circuit 43 on its negative input. 
30 Hence, the subtractor 45 produces a value equal to the difference between the sum of the squares 
of the samples and nji . This corresponds to the parenthetical portion of expression (10). 



14 



[0074] The value from the subtractor 45 goes to one input of yet another multiplier 

circuit 47 or other two-input product device. The multiplier circuit 47 receives the value for 
1/(1 -n) from the circuit 37. Hence, the multiplier circuit 47 computes the value of the variance 
(a 2 ) in accord with the expression (10), for output from the second output lead (Out2) of the 
module 19. An arithmetic logic unit performs the math function 49 for computing the square 
root of the variance and supplies the resultant value for the standard deviation (a) on the third 
output lead (Out3) of the module 19. 

[0075] Those skilled in the art will recognize that the present concepts have a broad 

range of applications, and the embodiments admit of a wide range of modifications, without 
departure from the inventive concepts. For example, the illustrated embodiment provided the 
mean, the standard deviation and the variance. It is a simple matter to add further computation 
modules to process one or more of these statistical values to produce desired measures of 
performance, such as Cpk. 

[0076] The embodiments may utilize variable samples taken automatically during 

operation of a machine, for example a disk drive or tape drive. Alternatively, the inventive 
technique may utilize variable samples taken by testing or monitoring during a manufacturing 
process, for example, by sampling the performance of products created by the manufacturing 
process. In either case, the inventive technique accumulates the necessary sums in real-time 
during the actual sampling, and without longer term storage of the original sample data. When 
all samples have been taken, a simple set of one-time calculations provides the mean, the 
standard deviation and/or the variance. More complex processing of stored individual samples 
off-line, after completion of the sampling routine, is no longer necessary. 

[0077] One advantage of the inventive real-time technique is the ability to imbed the 

calculation into an electronic product, itself For example, the computation capability may be 
imbedded in a tape or disk drive. When instructed to perform a test, for example during 
certification testing immediately following manufacture, the processor of the drive will calculate 
and output the desired values in real-time, substantially at the end of the time needed to measure 
the desired number of samples of the test parameter. The data need not be processed off-line by 
a separate computer. 
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[0078] Fig. 4 shows the functional elements of a data device, such as a disk drive or a 

tape drive. In this example, the device or drive 50 has the exemplary embodiment of Figs. 1-3 
embedded therein. 

[0079] The drive 50 includes a variety of physical elements 51 for handling the media, 

such as a tape or disk, and for reading and/or writing data to the media. Among or associated 
with the drive elements 51, the drive includes one or more sensors 53, for measuring at least one 
significant parameter during operation of the drive elements. The actual measurement signal 
from the sensor 53, in this example, is assumed to be an analog signal, although in some devices 
the signal may be a digital signal. If the sensor 53 provides an analog signal, that signal goes to 
an analog-to-digital converter (ADC) 55. 

[0080] The ADC 55 sequentially supplies samples of the measured signal to a 

micrcontroller 57, at a rate determined by the sampling clock 59. The clock 59 also provides 
timing signals to the micrcontroller 57. Although shown separately for discussion purposes, the 
clock 59 may be incorporated into the micrcontroller 57. 

[0081] The drive also includes a random access memory (RAM) 61 or other dynamic 

storage device, coupled to the microcontroller 57 for storing information as used and processed 
by microcontroller 57. The RAM memory 63 also may be used for temporary storage of 
executable program instructions. The drive 50 further includes a program memory 63, for 
storing the microcode program for the microcontroller 57. The memory 63 typically comprises 
read only memory (ROM) and/or electrically erasable read only memory (EEROM). In accord 
with the inventive embodiments, the microcode stored in the memory 63 includes the 
programming for the interim computation module 17 and the programming for the one time 
computation module 19. 

[0082] In a known manner, the microcontroller 57, itself comprises registers and other 

components for implementing a central processing unit and possibly an associated arithmetic 
logic unit. In operation, the accumulations performed in the module 17 may utilize registers in 
the control processing unit or appropriate sections of RAM memory 63. The microcontroller 57 
may be a microprocessor, a digital signal processor or other programmable device, implemented 
either as a general purpose device or as an application specific integrated circuit (ASIC) chip. 
Although shown separately, some or all elements of the memories 61, 63 may be incorporated 
directly within micrcontroller 57. 



[0083] The central processing unit of the microcontroller 57 executes the microcode for 

the modules 17 and 19 as well as the sample counter 21 and the constant value (n) input 23, to 
perform the computations of mean (/a), variance (a 2 ) and standard deviation (a), exactly as 
described above relative to Figs. 1-3. The microcontroller 57 outputs these statistical measures 
of the performance of the drive 50 via a data port 65, which may be the same port through which 
it outputs and receives normal data going to or from the medium. Alternatively, the 
microcontroller 57 may output the statistical measures of drive performance via one or more 
separate test data ports. The measures of performance typically go to a test computer or the like, 
for display and/or further processing. If the drive is already installed at the time of the test, the 
various results would go to the host computer for display and/or further processing 
[0084] As shown by the discussion above, the embodiments may be implemented using 

discrete logic circuits. Preferably, embodiments may be formed by appropriately programming a 
processor or microcontroller to process the measured parameter samples as taught above. Hence, 
some embodiments of invention encompass programmed processors for implementing the 
statistical computations in the inventive manner. Other embodiments of invention encompass a 
software product in the form of at least one machine-readable medium and one or more modules 
of code carried by the medium that are executable by a processor to cause the processor to 
compute one or more statistical measures of performance in essentially the manner performed in 
the embodiments discussed and illustrated herein. 

[0085] A machine-readable medium, as used herein, may be any physical element or 

carrier wave, which can bear instructions or code for performing a sequence of steps in a 
machine-readable form or associated data. Examples of physical forms of such media include 
floppy disks, flexible disks, hard disks, magnetic tape, any other magnetic medium, a CD-ROM, 
any other optical medium, a RAM, a ROM, a PROM, an EPROM, a FLASH-EPROM, any other 
memory chip or cartridge, as well as media bearing the software in a scannable format. A carrier 
wave type of medium is any type of signal that may carry digital information representative of 
the data or the instructions or code for performing the sequence of steps. Such a carrier wave 
may be received via a wireline or fiber-optic network, via a modem, or as a radio-frequency or 
infrared signal, or any other type of signal which a computer or the like may receive and decode. 
[0086] While the foregoing has described what are considered to be the best mode and/or 

other preferred embodiments, it is understood that various modifications may be made therein 



and that the invention or inventions disclosed herein may be implemented in various forms and 
embodiments, and that they may be applied in numerous applications, only some of which have 
been described herein. It is intended by the following claims to claim any and all modifications 
and variations that fall within the true scope of the inventive concepts. 



